Montonic convergent optimal control theory to modulate bandwidth limited laser 
pulses in linear and non-linear optical processes 
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We present a modified optimal control scheme based on the Krotov method, which allows for strict 
limitations on the spectrum of the optimized laser fields, without losing monotonic convergence of 
the algorithm. The method guarantees a close link to learning loop control experiments and is 
demonstrated for the challenging control of non-resonant Raman transitions, which are used to 
implement a set of global quantum gates for molecular vibrational qubits. 
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With the progress of laser pulse shaping and learning 
loop techniques [l|| quantum control experiments (OCE) 
became a forefront tool for the control and deciphering of 
molecular quantum processes 0,i,S,li- Optimal con- 
trol theory (OCT) [q, 0| as the theoretical counterpart 
is a powerful method for the prediction of pulse struc- 
tures as initial guess and guidance for OCE. With OCT, 
insight into the quantum pathways of these processes is 
directly available. The numerous applications of optimal 
control range from the control of chemical reactions in 
gas and condensed phase (8|, [oj to the control in nanos- 
tructures (lol . [ll| and to quantur a op tical problems like 

or the 



quantum information processing 



12, 
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preparation of cold molecules 3, 

One fundamental diflFerence between OCE and OCT 
is the spectral bandwidth of the laser field inherently 
present in the experiment but in principle unlimited in 
the original theoretical formulation. The general compa- 
rability of experimental and theoretical results may be 
complicated, since the theoretical answer for the optimal 
pulse can always span a wide bandwidth with quantum 
pathways out of experimental reach. Several suggestions 



suggestions 



have been made dealing with this challenge 
however, at the cost of monotonic convergence or general 
applicability [21]. 

We present a modified OCT approach based on the 
Krotov method [22] that treats time and frequency do- 



main equally while providing monotonic convergence. 
The method oflFers an elegant possibility to study OCEs 
theoretically by explicitly including as a constraint the 
crucial experimental feature of the spectral bandwidth. 

As an ultimate test we demonstrate the new tool by 
the implementation of stimulated non-resonant Raman 
quantum gates for vibrational qubits. The idea of molec- 
ular vibrational quantum computing [l2| has first been 
introduced for IR- active modes. Ultrashort, specially 
shaped laser pulses act as global quantum gate oper- 
ations. DiflFerent types of quantum gates and quan- 
tum algorithms have been demonstrated theoretically for 
IR transitions [H, [13, S [H, ^ and STIRAP pro- 
cesses (27|,l28|. Experimentally, molecular quantum gates 
have been realized in the^sible regime fiol, l30l| and latest 



IR shaping experiments [3l|, l32| open the route for the re- 
alization of vibrational qubits in the IR [33]. Stimulated 
non-resonant Raman quantum gates will provide new 
fiexibilities, like the choice of laser wavelengths. Their 
theoretical implementation can be regarded as a great 
challenge for the new OCT scheme since it comprises a 
non-resonant, two-photon, two-color process. 

The multi-target optimal control functional [ill for the 
molecular non-resonant Raman interaction includes two 
laser fields ei{t) with / = 1,2 (Eq. [T]). The corresponding 
Raman Hamiltonian is given by Vr = — ^ei(t) a e2(t), 
with the molecular polarizability a. 
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The objective is determined by the square of each over- 
lap (^/e(T)|$/e) of the propagated initial states ^fe(^) at 
the final time T with the target states ^/c of the global 



quantum gate operation. The change of the pulse energy 
is restricted with the Krotov-change parameter a^. ei(t) 
are the reference fields. A temporal shape function |[3^] 
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denoted as s{t) is inserted to achieve smooth switching 
on and off behavior of each laser field. The wave function 
has to satisfy the time dependent Schrodinger equation 
including the time evolution of the non-resonant Raman 
process. With the new term Fi{ei{t)) we introduce a fre- 
quency filter operation in its time representation, which 
restricts each electric field ei{t). The corresponding La- 
grange multipliers are 7z(t). The filter operations can in 
principle be realized in the time domain by linear digital 
filters, and particularly by finite impulse response (FIR) 
filters H 



N 



F{e{t)) = ^Cje{t-jAt), 



with the FIR filter coefficients cj and the step size At 
in the discrete time representation. By variation of the 
functional (Eq.[T]) with respect to the initial states ^/c(t), 
the target states ^k{t) and the laser fields ei{t) a set of 
coupled differential equations can be derived. The itera- 
tive calculation of the laser fields is performed with the 
Krotov method [22] and the next iteration step n + 1 
for the laser field ei{t) and analogously for 62 (t) can be 
formulated as: 



and 



with e1(t) = ii(t) 



(2) 



(3) 



The Lagrange multipliers 7/ {t) can be interpreted as cor- 
rection fields needed to suppress the undesired frequency 
components. In the optimal case the Lagrange multipli- 
ers ji{t) are adjusted to substract exactly the undesired 
field components from the optimized uncorrected fields 
^j^Ci^k (Eq. [3]). The spectral constraint Fi{ei{t)) de- 
pends only linearly on each electric field and it is possible 
to realize the side conditions \Fi{ei{t))\ = using Fourier 
filters fi{u;). It turned out that under practical consider- 
ations it is easier to use Fourier filters instead of FIR fil- 
ters. The Lagrange side conditions can be implemented 
in form of band-stop filter operations using the inverse 
//(cj) = 1 — fi{u;) of the band-pass filters fi{uj) which 
guarantees that only the undesired spectral components 
pass the band-stop filter. The Lagrange multipliers 7^(t) 
cannot be determined directly. In fact, for the calculation 
of 7/(t) the field change Ci^k must be predicted in the 
actual iteration step. This task is performed by propa- 
gating the target states and the intital wave functions 
with the laser fields from the previous iteration. 
The construction of the resulting fields j^t) resembles 



the OCT fields of the unmodified algorithm. 
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Filtering this output 7i (t) or analogously 72 (t) with the 
band-stop filter /[(oo) transforms them into the correction 
fields 71 (t) or 72 (^). The transformation is accomplished 
with the help of Fourier transforms JT. 



7jW = .^-M//H-WW)] 



(4) 



Each of the new fields e[^^^(t) can now be calculated by 
inserting the result from Eq. [H in Eq. O The Lagrange 
multipliers 7^(t) represent the time dependent electric 
fields of the undesired frequency components. The cor- 
rection fields are evaluated in each iteration step and are 
substracted from the optimized uncorrected field (Eq.[2j). 
Finally, to maintain the validity of the side condition the 
optimized field has to be filtered with the band-pass fi{uj) 
after each iteration. 

The modified OCT-scheme provides monotonic conver- 
gence, i.e. each iteration step improves the objective. Its 
convergence is proved in close analogy to the procedure 
given in fl6'| for the standard Krotov OCT. The differ- 
ence in the line of argumentation arises from the new 
constraints and enters in: 
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dt > 0, 
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corresponding to Eq. A15 of [l^. Aei{t) 
ef{t) denotes the change of the laser fields between two 
iterations. 

In accordance with a Lagrange side condition the out- 
put of the filter operations has to be zero. Consequently, 
all terms of Eq. [5] containing Fi{e{t)) vanish. Inserting 
Eq. Elinto Eq. [1] leads to: 



^3s{t) 
4 ao 




dt > 0, 



where only one non-quadratic term appears besides two 
positive, quadratic ones. As a result, the integral is al- 
ways greater than or equal to zero, which meets the re- 
quirements of convergency. We demonstrate the strength 
of the proposed OCT-algorithm with the implementa- 
tion of highly efficient stimulated non-resonant Raman 
quantum gates. A schematic sketch of the vibrational 
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FIG. 1: Stimulated non-resonant Raman quantum gates with 
the four two-qubit basis states, a) The exact manipulation of 
all quantum transitions indicated by the arrows (|00) ^ |01) 
and 1 10) ^ 1 11)) corresponds to a global NOT gate operation. 
A CNOT gate is realized by pulses only switching the state 
of the second (active) qubit when the first (control) qubit is 
in state |1). b) OCT optimization of a single Raman field 
without frequency restrictions leads to a spectrum with an 
additional carrier frequency 002 separated by 2A with respect 
to 002. 



ladder and the two quantum gate operations, NOT and 
controlled-NOT (CNOT), is depicted in Fig.[T](a). Stan- 
dard OCT-schemes fail for the simultaneous optimiza- 
tion of two non-resonant laser fields, since the virtual 
states are not determined within the formalism and the 
carrier frequencies of the laser pulses are independent 
of the eigenvalues of the system. Moreover, the elec- 
tric fields ei{t) and 62 (t) are not distinguishable. When 
one laser field of the Raman process is kept constant, 
the frequency unrestricted optimization of the other field 
yields a second band UJ2 in its spectrum (Fig.[T](b)). The 
additional band is shifted to higher frequencies by two 
times the energy difference A of the addressed qubit basis 
states. Such a broad spectrum requires extremely short 
laser pulses. The new algorithm provides the opportunity 
to optimize both laser fields within a selected and lim- 
ited frequency range and simple structured, stimulated 
non-resonant Raman quantum gates with high efficien- 
cies could be predicted for the first time. 

As a model system we selected two strongly Ra- 
man active C-H stretching vibrations of n-butylamine. 
The potential energy surface as well as the polariz- 
ability tensor components were calculated with density 
functional theory (b31yp/6-31++G**) ^] along both 
modes. The eigenfunctions and eigenvalues were ex- 
plicitly evaluated by a relaxation method (Hi. The 
quantum dynamics were carried out with a Cheby- 
chev propagation scheme (33|. Both selected modes 
with the fundamental frequencies ui = 2990 cm ~^ and 
U2 = 3030 cm ~^ provide high but balanced anharmonic- 
ities (intramode Ai = 74 cm A2 = 103 cm ~^ and in- 
termode Ai2 = 22cm~^), which are favorable molecu- 
lar properties for vibrational quantum computing | >23i] . 
For the definition of the two-qubit basis (|00), |01), |10), 
|11)) as sketched in Fig. [T](a) we encode the vibrational 
ground state of each selected normal mode as the logic 
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FIG. 2: The Raman fields ei(t) are depicted, 62 (t) have the 
same envelope functions, but different carrier frequencies, a) 
Laser field for the global CNOT. b) Spectra of both CNOT 
fields with their respective band-pass functions (dashed lines), 
c) Laser field for the global NOT. 



value and the first excited state as the logic value 1. 
The OCT-calculations were performed in the eigenstate 
representation, using the 50 lowest eigenstates. For the 
description of the laser-molecule-interaction, we selected 
the x^-tensor component surface and evaluated the cor- 
responding matrix elements. A universal set of quan- 
tum gates is implemented for this two-qubit system by 
stimulated non-resonant Raman processes. The CNOT 
and NOT gate with efficiencies above 99 % are exemplar- 
ily presented. ei{t) and 62 (t) of the global CNOT gate 
can be realized by simple gaussian-shaped laser fields 
(compare Fig. [2](a) for ei{t)). Their related spectra 
are depicted together with their band-pass filter func- 
tions fi{u;) in Fig.[2](b). The carrier frequencies were cho- 
sen to be in the near IR-regime with 800 nm (12500 cm~^) 
and 643 nm (15541 cm-^). The global NOT operation 
(Fig. [2](c)) is more complex in structure, because two 
vibrational transitions (compare Fig. [T](a)) have to be 
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FIG. 3: The solid line indicates the convergence of the mod- 
ified OCT algorithm after the total time T (determined by 
the yield Xlt=i during the optimization. 

The amount of the suppressed frequency components during 
the optimization is shown by the dashed line (normalized to 
unity). 
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driven simultaneously and the fundamental qubit tran- 
sition is close to the corresponding passive one. The 
x^-tensor component of the polarizability drives the vi- 
brational transitions for both qubit modes. Thus, it can 
be expected that polarized fields shaped to discriminate 
the qubit modes, might further simplify the laser field 
structures. 

The monotonic convergence for the CNOT gate op- 
timization with a Krotov-change parameter ao=10 can 
be traced from Fig. [31 solid line. The dashed line shows 
the evolution (normalized to unity) of the undesired spec- 
tral components during the optimization. Since the guess 
fields were chosen as simple, bandwidth tailored gaussian 
shaped laser fields, no frequency components have to be 
suppressed. In the first iteration step the amount of un- 
desired spectral components jumps to a maximum of less 
than 1 % of the pulse energy, but converges to zero while 
reaching the optimization aim. 

In summary, we present a Krotov OCT approach, 
which treats time and frequency domain equally, thus, 
unifying global optimal control with spectral constraints. 
The new tool optimizes laser fields under realistic ex- 
perimental spectral conditions. Optimal laser fields and 
control pathways in the experimentally accessible search 
space will be predicted. Additionally, an arbitrary pat- 
tern can be imprinted on the selected frequency range to 
suppress or enhance distinct quantum pathways. Thus, 
a strong and direct link to OCE is provided. The 
method has been successfully demonstrated for a multi- 
photon process and can easily be transferred to linear 
processes. A new realization strategy for molecular vibra- 
tional quantum computing has been presented by the im- 
plementation of simple structured, non-resonant stimu- 
lated Raman quantum gates of high efficiencies. Thereby, 
vibrational quantum computing in the ground state is 
transferred for the first time to the 800 nm regime, which 
is well established and accessible for pulse shaping tech- 
niques. 

The Deutsche Forschungsgemeinschaft through the ex- 
cellence cluster Munich Centre for Advanced Photonics 
and through the Normal verfahren is acknowledged. 



[7] W. Zhu, J. Botina, H. Rabitz, J. Chem. Phys. 108, 1953 
(1998). 



[1] P. Nuernberger, G. Vogt, T. Brixner, G. Gerber, Phys. 

Chem. Chem. Phys, 9, 2470 (2007). 
[2] A. Assion, T. Baumert, M. Bergt, T. Brixner, B. Kiefer, 

V. Seyfried, M. Strehle, G. Gerber, Science 282, 919 

(1998). 

[3] H. Rabitz, R. de Vivie-Riedle, M. Motzkus, K. Kompa, 

Science 288, 824 (2000). 
[4] J.L. Herek, W. Wohlleben, R.J. Cogdell, D. Zeidler, M. 

Motzkus, Nature 417, 533 (2002). 
[5] C. Daniel, J. Full, L. Gonzalez, C. Lupulescu, J. Manz, 

A. MerU, S. Vajda, L. Woste, Science 299, 536 (2003). 
[6] A. Bartana, R. Kosloff, D.J. Tannor, J. Chem. Phys. 106, 

1435 (1997). 



T. Brixner, G. Gerber, Chem. Phys. Chem. 4, 418 (2003). 
R.J. Levis, H.A. Rabitz, J. Phys. Chem. A 106, 6427 
(2002). 

M. Aeschlimann, M. Bauer, D. Bayer, T. Brixner, F.J. 
Garcia de Abajo, W. Pfeiffer, M. Rohmer, C. Spindler, 
F. Steeb, Nature 446, 301 (2007). 

J. Giidde, M. Rohleder, T. Meier, S.W. Koch, U. Hofer, 
Science 318, 1287 (2007). 

CM. Tesch, R. de Vivie-Riedle, Phys. Rev. Lett. 89, 
157901 (2002). 

R. de Vivie-Riedle, U. Troppmann, Chem. Rev. 107, 
5082 (2007). 

A. Sporl, T. Schulte-Herbriiggen, S.J. Glaser, V. 
Bergholm, M.J. Storcz, J. Ferber, F.K. Wilhelm, Phys. 
Rev. A 75, 012302 (2007). 

U. Dorner, T. Calarco, P. Zoller, A. Browaeys, P. Grang- 
ier, J. Opt. B: Quantum Semiclass. Opt. 7, S341 (2005). 
CP Koch, J.P Palao, R. Kosloff, F. Masnou-Seeuws, 
Phys. Rev. A 70, 013402 (2004). 

D.J. Tannor, A. Bartana, J. Phys. Chem. A 103, 10359 
(1999). 

M. Artamonov, T.-S. Ho, H. Rabitz, Chem. Phys. 305, 
213 (2004). 

P. Gross, D. Neuhauser, H. Rabitz, J. Chem. Phys. 96, 
2834 (1992). 

J. Werschnik, E.K.U. Gross, J. Opt. B: Quantum Semi- 
class. Opt. 7, S300 (2005). 

T. Hornung, M. Motzkus, R. de Vivie-Riedle, J. Chem. 
Phys. 115, 3105 (2001). 

J. Somloi, V.A. Kazakov, D.J. Tannor, Chem. Phys. 172, 
85 (1993). 

C Gollub, U. Troppmann, R. de Vivie-Riedle, New J. 
Phys. 8, 48 (2006). 

U. Troppmann, C Gollub, R. de Vivie-Riedle, New J. 
Phys. 8, 100 (2006). 

T. Cheng, A. Brown, J. Chem. Phys. 124, 034111 (2006). 

M. Ndong, D. Lauvergnat, X. Chapuisat, M. Desouter- 

Lecomte, J. Chem. Phys. 126, 244505 (2007). 

D. Sugny, M. Ndong, D. Lauvergnat, Y. Justum, M. 

Desouter-Lecomte, J. Photochem. Photobiol. A 190, 359 

(2007). 

C Menzel- Jones, M. Shapiro, Phys. Rev. A 75, 052308 
(2007). 

R. Zadoyan, D. Kohen, D.A. Lidar, V.A. Apkarian, 
Chem. Phys. 266, 323 (2001). 

J. Vala, Z. Amitay, B. Zhang, S.R. Leone, R. Kosloff, 
Phys. Rev. A 66, 062316 (2002). 

S-H. Shim, D.B. Strasfeld, M.T. Zanni, Optics Express 
14, 13120 (2006). 

D.B. Strasfeld, S-H. Shim, M.T. Zanni, Phys. Rev. Lett. 
99, 038102 (2007). 

B. M.R. Korff, U. Troppmann, R. de Vivie-Riedle, J. 
Chem. Phys. 123, 244509 (2005). 

K. Sundermann, R. de Vivie-Riedle, J. Chem. Phys. 110, 
1896 (1999). 

A.V. Oppenheimer, R.W. Schaffer, J.R. Buck, Discrete 
Time Signal Processing second editon. Prentice Hall (Up- 
per Saddle River, NJ 1999). 

M.J. Frisch et al., GAUSSIAN03, revision D.Ol, Gaus- 
sian, Inc., Wallingford, CT, 2004. 

T.H. Ezer, R. Kosloff, J. Chem. Phys. 81, 3967 (1984). 



